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Abstract: A new class of Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmological 
models with time-evolving fundamental parameters should emerge naturally from a de- 
scription of the expansion of the universe based on the first principles of quantum field 
(NJ ■ theory and string theory. Within this general paradigm, one expects that both the gravi- 

tational Newton's coupling G and the cosmological term A should not be strictly constant 
Tj- , but appear rather as smooth functions of the Hubble rate H(t). This scenario ("running 

\ FLRW model") predicts, in a natural way, the existence of dynamical dark energy without 



invoking the participation of extraneous scalar fields. In this paper, we perform a detailed 
study of some of these models in the light of the latest cosmological data, which serves 
to illustrate the phenomenological viability of the new dark energy paradigm as a seri- 
ous alternative to the traditional scalar field approaches. By performing a joint likelihood 



analysis of the recent supernovae type la data (SNIa), the Cosmic Microwave Background 
(CMB) shift parameter, and the Baryonic Acoustic Oscillations (BAOs) traced by the Sloan 
Digital Sky Survey (SDSS), we put tight constraints on the main cosmological parameters. 
Furthermore, we derive the theoretically predicted dark-matter halo mass function and the 
corresponding redshift distribution of cluster-size halos for the "running" models studied. 
Despite the fact that these models closely reproduce the standard ACDM Hubble expan- 
sion, their normalization of the perturbation's power-spectrum varies, imposing, in many 
cases, a significantly different cluster-size halo redshift distribution. This fact indicates that 
it should be relatively easy to distinguish between the "running" models and the ACDM 
using realistic future X-ray and Sunyaev-Zeldovich cluster surveys. 
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1. Introduction 

The high-quality observations performed during the last two decades have enabled cosmo- 
logists to gain substantial confidence that modern cosmology is capable of quantitatively 
reproducing the details of the many observed cosmic phenomena y], ||] , including the 
accelerated expansion of the universe at the present epoch The field of cosmology, 
therefore, is no longer a pure realm of philosophical speculation, but a rigorous branch 
of experimental physics. In view of this fact, the status of modern cosmology can be 
declared as very healthy and highly satisfactory from the point of view of the empirical 
sciences. Nevertheless, if we pause for a moment and make a deeper reflection, we soon 
come to realize that the ultimate understanding of these magnificent observations is still 
highly blurred, if not completely opaque, to the light of the first principles of fundamental 
theoretical physics, e.g. from the viewpoint of Quantum Field Theory (QFT) and string 
theory. 

Indeed, we realize that what we have been able to do throughout the last decade 
is to furnish a rather precise numerical fit to the parameters of the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) model in the light of the most accurate observations, but at the 
moment these observational successes are not accompanied by a deeper theoretical insight 
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into the meaning of the preciously collected data on dark matter (DM) and dark energy 
(DE). Needless to say, the resulting high quality of this combined fit (from independent 
data sources) speaks very well for the FLRW model as a likely correct starting point to 
understand the structure of our cosmos. However, it does not provide a single clue to the 
meaning of - or on the physical substratum lying underneath - the fitted energy density 
parameters at the present time and (except for the baryonic part — 0.04, 
of course, which we know from primordial nucleosynthesis), nor does it explain why the 
geometric parameter i s comparatively much smaller. In practice, observations tell us 
that Q°j ~ 0.28 and £l° A ~ 0.72 (with Q° K < 0.01. It follows that most of the matter content 
of the universe is non-baryonic DM, and that the value of the cosmological constant density 
p\ = A/(8nG), or, in general, the current dark energy (DE) density, is 

p\ = n° A pi ~ 6 h 2 x 10- 47 GeV 4 ~ 3 x 10" 47 GeV 4 , (1.1) 

where p° c ~ (s.O y/h x 10~ 12 GeV) ~ 4 x 10 47 GeV 4 is the present value of the critical 
density and Hq = 100 h km/s/Mpc, with h — 0.704, the current Hubble's rate. 

Despite the "scandal" of our having to admit full ignorance about the identity of most 
(~ 90%) of the matter budget of the universe, this trouble pales in comparison to the 
colossal enigma hidden behind the value (|Ll|). Understanding this value is perhaps the 
biggest scientific mystery of all times: "the cosmological constant problem" || |||, which 
manifests itself as the double conundrum of the tiny observed value of the cosmological 
constant (CC) in Einstein's equations - the "old CC problem" [0] -, and also the puzzling 
fact that this value is so close to the current matter density (the "cosmic coincidence 
problem" (8|), including the mystery of a possible phantom behavior near our time. Solving 
these problems surely requires introducing new fundamental ideas of profound theoretical 
meaning and scope. Some promising new avenues have recently been proposed ||, [l0|, [yj 
12, 13 1 1 - see also the recent alternative attempts [15, [[(], 17, 18 1. 



For a long time model builders have pursued an explanation of these cosmological puz- 
zles in terms of dynamical scalar fields [|l!j] , later adapted into the idea of quintessence . 
Quintessence was proposed without attempting to explain the smallness of the CC, but 
only to cope with some aspects of the cosmic coincidence problem ||. These scalar field 
models ignore ah initio the huge vacuum energy from the Standard Model (SM) of strong 
and electroweak interactions, which is 55 orders of magnitude larger than the measured 



value fll.lD — see |J (section 2 and Appendix B) for a very detailed discussion. Moreover, 
quintessence models usually postulate a preposterously tiny mass for these scalar fields of 
the order of 10~ 33 eV, i.e. comparable to the current Hubble rate Hq. This is a great 
disadvantage, not to mention the fact that they are plagued with the same fine tuning 
problems as the original CC approach itself, as emphasized in [j|]. 

From our point of view, a sensible possibility that should not be neglected is to think 
of a "dynamical A" or "effective vacuum energy" p\ e s, where, rather than replacing A by 
a collection of ad hoc ersatz fields, we stick to the idea that the CC term in Einstein's 
equations is still a "true cosmological term" , although we permit that "the observable CC 



1 For a recent review of some of these ideas, see |TJ|] and references therein. 
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at each epoch" can be an effective quantity evolving with the expansion of the universe: 
PAeff = PAes(H), where H is the Hubble rate. While the old phenomenological models 
for time-evolving CC [21, 22] lacked of a fundamental motivation, a general dynamical A 
approach within QFT has more recently been emphasized in [23, 24, 25, 26] and also in the 
past from the point of view of QFT in curved space-time - see |[^] and [27, 28, [29| ]. 

The more modern idea of a slowly running CC is based on the possibility that quantum 
effects in curved space-time pH] can be responsible for the renormalization group (RG) 



running of the vacuum energy. One can show that these models [ 14 [ are nicely compatible 
with the most recent experimental data [31], and this fact spurs us to further explore their 
impact both in the cosmological and astrophysical domain. We shall concentrate here on a 
class of RG models in which the vacuum energy is a function of cosmic time t though the 
Hubble rate H = H{t). An archetype example is the quantum field vacuum model based on 
the evolution law A(H) = uq + ni H 2 . In Ref. |H|] some of us have investigated the global 
dynamics of this cosmological model (together with various alternative A(t) models), in the 
light of the most recent cosmological data. However, in that paper we assumed that the 
vacuum decayed into matter, whereas here we wish to explore the complementary possibility 
that matter is conserved thanks to a logarithmically variable gravitational coupling G = 
G(H). This is very appealing because, if matter is conserved, then one cannot run into 
potential problems related with the fact that the vacuum energy could decay in part into 
baryons and photons, which would be of course undesirable. We shall call such FLRW-like 
model, with G = G(H) and A = A(H), the "running FLRW model" (hereafter A^CDM), 
in contrast to the "running ACDM model" (denoted A^CDM) in which only A evolves with 
time at the expense of maintaining an interaction with matter pi] . 

In this paper, we wish to test the two basic types of running models in the light of 
the most recent cosmological data on SNIa+BAO+CMB (i.e., distant type la supernovae, 
baryonic acoustic oscillations and cosmic microwave background anisotropies). The clus- 
tering properties of the vacuum energy can also be of high interest as they may help to 
shed some light on the fundamental issue of non-linear structure formation. The subject 



has been discussed for different DE models, although mainly based on scalar fields |32[]-[33] 
(and references therein). A first study in this direction for running vacuum models was 
recently presented by some of us in[^]. In the current work, however, we emphasize on 
the implications for the formation and distribution of the collapsed cosmic structures (i.e. 
galaxy clusters) along the lines of the methods utilized in |3l|, which we apply here to the 
A^GiCDM model, and leave the presentation of the clustering effects in this model for a 
separate work. Therefore, we focus here on the study of the dark-matter halo mass func- 
tion and the redshift distribution of cluster-size halos, and compare with the corresponding 
results for the AjG^CDM and the concordance ACDM model. 

The paper is organized as follows. In the next section, we discuss some scenarios 
with variable A and G. In Sect. |3|, the FLRW model with running G and A (i.e. the 
A^GfCDM model) is introduced, while in Sect. || we review the AfCDM model. The 
confrontation of these models with the latest cosmological data is made in Sect. |6|, and 
their predictions concerning the formation of the galaxy cluster-size halos and the evolution 
of their abundance is presented in Sect. f|. In the last section, we deliver our conclusions. 
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Finally, in an appendix we provide details of the calculation of the linearly extrapolated 
critical density for the RG models under consideration. 



2. Models with time-evolving cosmological parameters 

The idea that the vacuum energy should not be a 'rigid' quantity in cosmology is a most 



natural one [ 14 ] . It is difficult to conceive an expanding universe with a strictly constant 
CC energy density, p A = A./(8nG), namely one that has remained immutable since the 
origin of time. A smoothly-evolving vacuum energy p\ = p\(t) that inherits its time- 
dependence from cosmological functions such as the scale factor a(t) or the Hubble rate 
H(t), is not only a qualitatively more plausible and intuitive idea, but is also suggested 
by fundamental physics, in particular by quantum field theory (QFT) in curved space- 



time [3C]. To implement this notion, it is not strictly necessary to resort to ad hoc scalar 
fields, as usually done in the literature (e.g. in quintessence formulations and the like PQl). 
A "running" A term can be expected on very similar grounds as one expects (and observes) 
the running of couplings and masses with a physical energy scale in QFT. 

Therefore, the guiding paradigm that we adopt in this paper is that the CC density 
should naturally be a function of the cosmic time t, or equivalently of the cosmological red- 
shift z. The generalization of this possibility to all the cosmological parameters, with an eye 
on trying to explain many surprising features of the dark energy (DE) - such as the apparent 
phantom-like behavior near our time - has been put forward in |12[| . Within this context, all 
the cosmological parameters are viewed as effective couplings whose running can be studied 
from the semi-classical formulation of QFT in curved space-time |j(J. Interestingly enough, 
the possibility to derive this evolution from an effective action has also been addressed in 



some particular cases [23, 42]. Let us also remark that these ideas have been pursued in 



the literature from different points of view since long ago 43, 27, 25, 2£, [44[], and they have 
been recently re-examined both from the theoretical point of view [24] and also from the 
viewpoint of various kinds of possible phenomenological implications []3l| , 45, 41, 4€, 47]. 

In this framework, the variation of "fundamental constants" such as A and G could 
emerge as an effective description of some deeper dynamics associated to QFT in curved 
space-time, or quantum gravity or even string theory, all of which share the powerful renor- 
malization group (RG) approach. The latter entails the possibility to analyze the potential 
impact of the leading quantum effects from matter particles on the cosmological observables 
such as vacuum energy and gravitational coupling. These effects should eventually lead to 
definite time-evolution laws pa = Pi\{t) ,G = G(t). However, given a fundamental model 
based on QFT, the parameter p\ will primarily depend on some cosmological functions 
(matter density p m , Hubble expansion rate H, etc) which evolve with time or redshift. 
Similarly for the Newton coupling G. Therefore, in general we will have the two following 



functions of the cosmic time or the cosmological redshift |12[ : 

p A (z) = p A (p m (z),H(z),...), G(z)=G( Pm (z),H(z),...). (2.1) 

Of course, other fundamental parameters could also be variable. For example, the fine 
structure constant has long been speculated as being potentially variable with the cosmic 
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time (and therefore with the redshift) including some experimental evidences - see e.g. [48]. 
This might also have an interpretation in terms of the RG at the level of the cosmological 
evolution. However, in this paper we concentrate purely on the potential variability of 
the most genuinely fundamental gravitational parameters of Einstein's equations such as 
A and G. The possibility that G could be variable has been previously entertained in 
the literature, although most of these formulations are usually framed within the original 
Jordan and Brans-Dicke proposals p9| ] linked to the existence of a dynamical scalar field 
coupled to curvature. 

The cosmological constant contribution to the curvature of space-time is represented 
by the A g^ v term on the l.h.s. of Einstein's equations, which can be absorbed on the r.h.s. 
of these equations: 

R„u - \g^R = 8ttG r M „ , (2.2) 

where the modified energy- momentum tensor is given by T^ u = T^ IU + g^ lv pA- Here 
p\ = A/(8ttG) is the vacuum energy density associated to the presence of A, and T^ v 
is the ordinary energy-momentum tensor of isotropic matter and radiation. Modeling 
the expanding universe as a perfect fluid with velocity 4- vector field U^, we have T^ u — 
—Pm9^u + (Pm + PmjUfJJy, where p m is the proper isotropic density of matter-radiation 
and p m is the corresponding pressure. Clearly the modified defined above takes the 
same form as with p = p m + p\ and p = p m — pa- Therefore: 

Tf+u = ~P9fiu + (P + P)U^U U = (pA -Pm)g^u + {Pm+PmW^Uy . (2.3) 

With this generalized energy-momentum tensor, and in the spatially flat FLRW metric 
ds 2 = dt 2 — a 2 (t)dx 2 , the gravitational field equations boil down to Friedmann's equation 

and the dynamical field equation for the scale factor: 

a = — —G(p + 3p) a = — —G(p m + 3p m - 2 p A ) a . (2.5) 

Let us next contemplate the possibility that G = G(t) and p\ = pxif) can be both functions 
of the cosmic time. This is allowed by the Cosmological Principle embodied in the FLRW 
metric. The Bianchi identities (which insure the covariance of the theory) then imply 
(^GT^^j = 0. In the case of the FLRW metric, the previous identity amounts to the 
following generalized local conservation law: 

j t [G(pA+p m )] + 3GH(p m +p m ) = 0. (2.6) 



This equation is actually a first integral of the dynamical system ( j2.4J ) and ( |2.5| ), as can be 
easily checked. When G is constant, the identity above implies that p\ is also a constant, if 
and only if the ordinary energy- momentum tensor is individually conserved (xj^T^ LU = 0), 
i.e. p m + 3H (p m + Pm) = 0. This equation can be rewritten 

p'm(a) + ~ (Pm(a) + Pm(a)) = 0, (2.7) 
a 
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where the prime indicates differentiation with respect to the scale factor: /' = df /da for 
any function /(a). The solution of ( ^??| ) is 

Pm (a) = pi a -3(i-K*») = fP m (i + ,)3(i+^) (2 . 8) 

where is the current matter density and oj m = p m / ' p m is the equation of state (EoS) 
parameter for cold (w m = 0) and relativistic (uj m = 1/3) matter, respectively. We have 
expressed the result ( |2.8D in terms of the scale factor a = a(t) and the cosmological redshift 
z = (1 — a)/a. 

A first non-trivial situation appears when pa = p/\_(t) but G remains still constant. 
Then equation ([2j]) boils down to 

PA + Pm + 3H(p m +p m ) = , (2.9) 

This scenario shows that a time-variable pa cosmology may exist such that transfer of 
energy occurs from matter-radiation into vacuum energy, and vice versa (see section 4). 
However, let us remark that the time evolution of p\ is still possible in a framework where 
matter is strictly conserved, i.e. such that equation ( |2.7D is maintained. Such scenario 
is highly desirable and is perfectly consistent with the Bianchi identity at the expense of 
a time- varying gravitational coupling: G = G(t). Indeed, it is easy to see that equation 
(|2.7|) is compatible with ( |2.6[) provided that the following differential constraint between 
the variable G and the variable pa is satisfied: 

(pm + PA)dG + Gdp A = 0, (2.10) 



where p m in this equation is given by (2.8). To solve the model in the basic set of variables 
(H(z), p m (z),p(z),G(z), pa(z)), still another equation is needed. In the next section we 
discuss the possibility that the fifth missing equation is a QFT-inspired relation of the form 
PA = Pa{H{z)). 

3. A FLRW-like model with running G and A: the A t G t CDM model 

Consider the class of models in which the vacuum energy and the gravitational coupling 
evolve as power series of some energy scale p, such that the rates of change are given 
respectively by 

d PM ST a ,, 2k d f 1 \ R ,, 2k tort 

r fc=0,l,2,... r v fc=0,l,2,... 

These expansions can be considered as purely phenomenological ansatzs or, if we aim at a 
more fundamental description, as emerging from the quantum field RG running of the cos- 



mological parameters p7| , |28| , 46 1. Such running should ultimately reflect the dependence 
of the leading quantum effects on some physical cosmological quantity ^ associated with p, 
hence p\ = pa(0- The physical scale £ could typically be the Hubble rate H, or even the 
scale factor a, which in most of the cosmological past also maps out the evolution of the 
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energy densities with iJ[23|. Since in both cases £ = evolves with the cosmic time, 
the cosmological term inherits a time-dependence through its primary scale evolution with 
H(t) or a(t). In this context, the r.h.s. of the above equations defines the /3-functions for 
the running of pa an d G~ l in QFT in curved space-time. The setting p, — > H naturally 
points to the typical energy scale of the classical gravitational external field associated to 
the FLRW metric, i.e. the characteristic energy of the FLRW "gravitons" attached to 
the quantum matter loops contributing to the running of p\ and G~ l in a semi-classical 
description of gravity. Coefficients A 2 k , B 2 k in these formulas receive contributions from 



boson and fermion matter fields of different masses Mj, and the series (3T) become expan- 
sions in powers of the small quantities H/Mi 



28L 29]. The fact that only even powers of 
H are involved is dictated by the general covariance of the effective action |J, H 2 ■ Aft er 
integrating the formulas ( |3.1| ), the result may be looked upon as a specific implementation 
of the general equations ( |2.1| ), where the expansion rate H is singled out as the main cos- 
mological variable on which the functions p\ and G~ 1 depend and in terms of which are 
expanded. 



It is important to note that the expansions (3.1) are correlated by the Bianchi identity 



( 2.10 ) and therefore they must be consistent with it, order by order. Let us consider them 
more closely. We start by first rewriting the expression for the RG of the CC term in a 
more explicit manner, as follows: 



1 



dp\(p) __ 

din /it (47r) 5 



v + E 



D; 



2n 2 /i 2 + % 4 ), (3.2) 



where Mi are the masses of the particles contributing in the loops, and Bj, Cj, .. are dimen- 
sionless parameters. We have omitted the Mf terms in the expansion since, as explained 
in p8| , p9| ], such terms would trigger a too fast running of the CC. It follows that only the 
"soft-decoupling" terms of the form ~ Mf p? remain in practice. The integrated form of 
fO) reads 



Pa(H) = n + n 2 H 2 



(3.3) 



where we have set p = H and neglected the smaller higher order terms, which are of 
0(H 4 ). The additive constant no in (3.3) must be essentially given by the current value of 
the vacuum energy, p A , but not quite. Notice that if n 2 = 0, the CC is strictly constant and 
then no just coincides with p A . However, if the H 2 term is to play a significant role while 
preserving the approximate ACDM behavior, it should neither be negligible nor dominant 
as compared to no- Being the latter the leading term in the series expansion, it must still 
be of order n ~ p A ; thus the largest masses Mj associated to n 2 ~ ^2i^iMf should be 
in the ballpark of the mass Mx associated to a typical Grand Unified Theory (GUT) near 



the Planck scale Mp (see below). From the structure of the expansion (3.2), with fx = H 
and Mj ~ Mx, we reconfirm that no other even power H 2n (not even H 4 ) can contribute 
significantly on the r.h.s. of equation (|3.3|) at any stage of the cosmological history below 



2 In practice, if one tries to fit the data with a time dependent CC term which is linear in the expansion 
rate, i.e. of the form A oc H, the results deviate significantly from the standard ACDM predictions ^p| . 
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Mx < Mp. It is particularly convenient to rewrite the coefficients of ( |3.3[) as follows: 

n = p° A -^M 2 H 2 , n 2 = ^-M 2 , (3.4) 

with 

(3-5) 

Here, Hq is the value of the Hubble rate at present. A most important (dimensionless) 
parameter of the present framework is v, given in equation ( |3.5| ). It defines the main coef- 
ficient of the /3-function for the running of the vacuum energy For v = 0, the vacuum 
energy remains strictly constant, p A = p A , whereas for non- vanishing v the evolution law 
Q3.3j ) can now be written as 3 

p A (H) = p° A + ^-M 2 (H 2 -H 2 ). (3.6) 

Let us stress that, in the QFT framework, v is naturally expected to be non-vanishing owing 
to the quantum effects of matter particles. The coefficients Bi in (|3,5|) can be computed 



from the quantum loop contributions of fields with masses Mi, see e.g. [23] for a specific 
framework. It is customary to write v in the compact form 

I2vr My y ' 

with M an effective mass parameter of order of the average mass of the heavy particles 



of the underlying GUT, including their multiplicities 28]. From the expression (3J3), it 
is clear that the natural range of that parameter is |^| <C 1; and it can be both positive 
or negative, since a = ±1 depending on whether bosons or fermions dominate in the 
sum of loop contributions ( |3.5[ ). For instance, if GUT fields with masses Mj near Mp do 
contribute, then \v\ < I/(127r) ~ 2.6 x I(U 2 , but in general we expect it to be smaller 
because the usual GUT scales, say Mx ~ 10 16 GeV, are not that close to Mp ~ 10 19 
GeV. By counting heavy particle multiplicities in a typical GUT (i.e. the total number of 
particles with masses Mi ~ Mx), a natural estimate lies in the range v = ICU 5 — I0~ 3 P3| . 
Clearly, the vacuum energy density Q3.6|) is normalized to the present value, i.e. 



p A (H ) = p A ^^-n A H 2 Mp. (3.8) 

Recall that the observed numerical value of this expression is given by equation (|I . 1|) , 
corresponding to Q A ~ 0.7. We remark that the parametrization (|3.6| ) satisfies the afore- 
mentioned condition that no is the leading term and is of order p A , and at the same time 
the correction term is of order M 2 H 2 , with M ~ ^J~uMp a large mass, even if v is as small 
as, say, |^| ~ 10~ 3 or less. 



3 It is very interesting to notice that this quadratic evolution law for the vacuum energy with the expansion 
rate has also been suggested recently by alternative QFT methods, see |16|. 



-8- 



What about the running equation for G 1 in equation (|3.1| )? A more detailed form 
for it reads: 

x ' i i i 1 

If we take again /i = H, then, since we have H <C M, for all the relevant epochs of the 
universe - at and below the GUT scale -, this series should converge extremely fast, the 
first term being the dominant one and the remaining terms are inessential. The first term 
is of course a contribution to the running of G _1 which is directly driven by the heavy 
particle masses Mi of the GUT, and therefore can be parametrized as a coefficient times 
the Planck mass. Thus, we are essentially led to a RGE of the type 

" f ^] =2vM 2 P + (3.10) 



dlnfi \G 

where the coefficient in front of Mp has been fixed by the condition that this solution 
is compatible with the Bianchi identity ( |2.10 ) after we adopt the solution ( |3.6[) for the 
evolution of the vacuum energy and set p = H, too, for the running G (see below). We 
shall show below explicitly this consistency, but let us first notice that the solution of (|3.1C| ) 
satisfying the boundary condition G(fi = Hq) = Gq = 1/Mp is 

We see that g(H) depends on the parameter v and that, to first order, we have g(H) = 
l + 0(u), i.e. G{H) = Gq{\ + 0{v)). Thus, v plays also (in a very evident manner here) the 
role of the /3-function for the RG running of G. For example, for v > the gravitational 
coupling decreases with the characteristic cosmological energy H (hence G increases with 
the expansion); it follows that for v > the coupling G is asymptotically free (as the 
QCD coupling at high energies), whereas for v < the gravitational coupling behaves as 
the electromagnetic one, i.e. increases with the energy (equivalently, G decreases in this 
case with the expansion). In both cases the running of G is extremely slow, since it is 
logarithmic and with a very small coefficient \v\ <C 1. Let us note that if we would have 
omitted the leading ~ Mf term on the r.h.s. of equation ( |3.9| ) - as we have done with the 
Mf term on the r.h.s. of equation ( |3.2[ ) - we would have obtained a still slower running of 
the gravitational coupling: 

1 

l + UH 2 ~Hl)/Ml 

where £ is a dimensionless coefficient. In deriving equation ( |3.12 ) we have assumed that 
the r.h.s. of ( |3.9| ) is dominated by the next-to- leading term ~ H 2 (for /x = H) under 
the assumption that the ~ Mf terms cancel. Since H 2 /Mp is utterly negligible, this new 
running law for G is even less sensitive to time evolution for the present universe. In 
this case, the corresponding evolution of the CC following from the Bianchi identity ( 2.1C| ) 
would also be very mild: 

p K {H) = pi + J- (H*-H*) . (3.13) 



9(H) = rr2WB<f2 , (3-12) 
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However, since there is no a priori theoretical reason nor a phenomenological constraint to 
assume that the leading term ~ Mf terms cancel, we shall concentrate here only on the 
form ( |3.11| ) and corresponding ( |3.6[ ), rather than on ( |3.12] )-( |3.13| ), when we confront our 
running model with the observational data in section H 

Before checking the consistency of the solution ( |3.11| ), it is convenient to introduce 
some notation. We define the energy densities normalized with respect to the current 
critical density = 3Hq / (8n Go): 

n^^PM- (i = m,A). (3.14) 

Pc 

It is also convenient to introduce the energy densities normalized with respect to the critical 
density at an arbitrary redshift, p c {z) = 3H 2 (z)/(8tt G(z)), in which both H and G are 
functions of z. This new set of normalized energy densities is given by 



Pc(z) E?(z) 



fii(z) (i = m,A), (3.15) 



where we have defined the normalized Hubble rate of the running FLRW model with respect 
to the current one: 

E(z) = ^ = VsR [M*) + ^a(^)] 1/2 , g(z) = ^ . (3.16) 

This is the generalized Friedmann's equation for the present model. Obviously, the pa- 
rameters ( 3.15 ) coincide with the previous ones (|3.14 ) only at z = 0, and at that point 



they all furnish the normalized current densities: f2j(0) = &i(0) = Q®. However, only the 
tilded cosmological parameters ( 3.15| ) satisfy the (flat space) cosmic sum rule (valid at any 



redshift z): 

Q m {z) + Q A {z) = 1 . (3.17) 

We agreed to call the FLRW-like model with variable G = G(H) and A = A(H) the 
A^GiCDM model, in contrast to the A^CDM model (previously studied in p8| and most 
recently in great detail in |pl[ ), in which only A evolves with H at the expense of main- 
taining an interaction with matter through the generalized conservation law ([?!]) - see the 
next section for a summarized discussion of the A^CDM model. 

The solution of the running FLRW model is determined from the Friedmann equation 
(3.16), the Bianchi identity ( |2.10| ) and the RG law for the CC ([O]). Using the normalized 



density parameters defined in equation ( 3.14 ), the basic cosmological equations can be 
formulated in the following compact way: 

E\z) = g(z) [n m (z) + n A (z)] , (3.18) 
(O m + n A )dg + g dQ A = , (3.19) 
n A (z) =n° A + v [E 2 (z) - 1] , (3.20) 

n m (z) = n m (i + z)^ 1+ ^\ (3.2i) 

where the last equation just reflects the covariant conservation of matter (for cold dark 
matter we have oj m = and for relativistic matter we use oo m = 1/3), i.e. it is a rephrasing 
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of equation (|2.8|) . Solving the differential form ( 3.19| ) in combination with the rest of the 
equations, it is easy to arrive at 

dq <Kl\ dE 

-4 = ± = -2v — . 3.22 

g 2 E 2 E V ' 

Integrating this expression with the boundary condition g(z = 0) = 1, and using also the 
fact that E(z = 0) = 1, we finally meet the explicit form ( |3.11| ) for the function g = g(H), 
in full consistency with our expectations — q.e.d.. 

We may also check a posteriori the consistency of the scale choice [i = H for the 
running of the cosmological parameters. The election of this scale is of course a difficult 



issue in cosmology, as there are different possibilities at our disposal - see [^, 23] and 
references therein. However, we expect that the choice [x = H should naturally be self- 
consistent in the FLRW context at the leading order. We can substantiate this assertion 
starting from the generic functions (|2.1| ), which we assume that are dominated by a single 
scaling variable fi, i.e. p\ = pa(m) an d G = G{(j). The differentials of these variables in 
the Bianchi identity fl2.10|) can then be written in terms of dfj,. We easily arrive at the 
following expression: 



dpA ( d 1 \ 1 3 2 



GM(p,„ + ft ,(M)) = ^^^j = s ^. (3.23) 
where in order to evaluate the final form on its r.h.s. we have used the running laws 



( |3.2| ) and (3.10) as well as the explicit form of the coefficient ri2 in fl3.4|) . Notice that 
the small parameter v has neatly canceled on the r.h.s. of ( |3.23| ), within the order under 
consideration. At the end of the day, equation ( 3.23j ) tells us that \x is determined by the 



simple formula 

M 2 = ^|^K + PA(M)], (3.24) 

which patently shows, upon comparison with (|2.4|), that \x = H for the flat space — q.e.d. 
Therefore, the identification [i = H is consistent within the order under consideration. 
Furthermore, it also suggests that the evolution laws ( |3.6| ) and ( p. 11 ) can be interpreted 



within the context of the RG in cosmology, with fi = H acting as the natural running scale. 

Let us now come back to the running of the gravitational coupling with fi = H. 
Although the form ( |3.11 ) is the simplest and most revealing one from the physical point 



of view, it is also convenient to determine the dependence of the gravitational coupling as 
a function of the cosmological redshift, i.e. g = g{z). This function is necessary for the 
numerical analysis of the model and the comparison with the observational data. To this 
end, we start by combining equations ( |3.18| ) and ( [3.20[ ), and we obtain: 

1 - v g{z) 



From the above equations we also have 



n m (z) + n A (z)= ' — . (3.26) 

1 - ug(z) 
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On comparing it with the sum rule ( 3,17| ) for the tilded parameters, we see that only for 



z = the r.h.s. of equation ( |3.26[ ) gives 1, as expected. Next we compute dftA(z) explicitly 
from equation ( |3.25| ) and substitute the obtained result, together with equation (|3.26 ), in 
the differential form ( |3.19| ). After some straightforward algebra, we find 

(n m (z) + fl A -v)dg + v{l-vg)g 2 dn m (z) = , (3.27) 

which can be integrated by quadrature: 

1 



' 1 + v In 



In [n m (z) + n° A - v] . (3.28) 



Notice that Q, m (z) in the previous equation is given by (^2l|) , but in those cases when 
there is a comparable mixture of cold matter and radiation we have to write 

n m (z) = n° m (i + zf + n R (i + z)\ (3.29) 

where = (1 + 0.227^) f2^, with N v the number of neutrino species and h? ~ 
2.47 x 10" 5 . This is the case, for example, when we analyze the matter content near the 
decoupling time for the analysis of the CMB observables, see the next section. 

Equation ( |3.28j ) determines g = g(z) as an implicit function of the redshift. There 
is no obvious way to write the explicit function. It is nevertheless easy to check that it 
satisfies the boundary condition g(0) = 1, as it should. With g{z) determined from fl3.2£j ). 
the redshift evolution of the vacuum energy follows from (|3.25| ) and ( |3.29| ) . For the present 
model, however, it is not possible to find an analytical expression of the cosmological 
functions with respect to the cosmic time, only in terms of the redshift, or equivalently the 
scale factor. But this is enough for the numerical analysis. 

As we know, the coefficient v in equation (|3.6| ), or equivalently ( |3.20| ), measures the 
amount of running of the CC or vacuum energy in the running FLRW model. For any 
given v, we can compare the value of &a(z) with the current value Q A ~ 0.7. The relative 
correction can be conveniently expressed as follows: 

AO A (.) - ^ A( l"^ = tSt [^(*) ~ 1] • (3-30) 

Since g = 1 for v = 0, it follows that, for small v, g(z) deviates little from 1, namely 
g(z) = 1 + 0(v). Thus, expanding to order v in the matter epoch, it is easy to show from 
the previous equation that 

O 

AOa(^^ [(l + ^) 3 -l] , (3.31) 

where g(z) ~ 1 to this order. If we look back to relatively recent past epochs, e.g. exploring 
redshifts z = 0(1) relevant for Type la supernovae measurements, the deviation AQa(z) 
is of order of a few times v. For example, for z = 1.5 and z = 2, the deviations are 
Ar?A(l-5) ~ 6f and AQ\(2) ~ 11 1/ respectively, assuming tt° m = 0.3. Although the 
correction is small, it is not necessarily negligible, even if \v\ ~ 10~ 3 . In the next sections, 
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we will see that the analysis of the cosmological observables in the light of the latest 
SNIa+BAO+CMB data amounts to a bound on \v\ of this order, and therefore there is 
hope to eventually becoming sensitive to the effects of the running vacuum energy. 

Since we will also use the data on structure formation in order to constraint the running 
A^GiCDM model in section |(| let us recall that in any model with variable G and p\ in 
which matter is covariantly conserved, the matter density contrast D = 5p m / p m satisfies 
the following second order differential equation with respect to the scale factor [EBfl: 



where the cosmological parameter O m (a) has been defined in ( 3.15| ). Moreover, the per- 



turbations in G are actually tightly linked to those in p\, and they are at the same time 
linked to the matter perturbations, as follows: 

5pA_JG 5 Pm _ (5G(a)Y 

PA ~ G ' Pm ~ (7(a) • ^ 

If we neglect the perturbations in G (5G ~ 0) and hence in p\ too, equation ( |3.32| ) reduces 
to the standard form of the growth factor for models with self-conserved matter (under the 
assumption of negligible DE perturbations), as for example in the case of matter pertur- 
bations in the ACDM model |33]]. If, on the contrary, we consider non- vanishing SG, and 
hence non- vanishing 5p\, one arrives at the following third order differential equation: 

D"'(a) + \ (l6 - 9^ m (a)) + | (s - im m (a) + (a) - a fi' m (a)) ^ = . 

(3.34) 

This equation was derived analytically (in both the synchronous and Newtonian gauges) 
in [45], where it was also numerically analyzed. Let us stress that when G is variable, 



it is necessary to include the perturbations in this variable since, as we have seen, they 
become correlated with the perturbations in pa because of the Bianchi identity of the 
running FLRW model, equation ( 2.10| ). Moreover, there is no consistent way to set them 



to zero if the matter perturbations themselves are nonzero, see equation ( |3.33| ). Therefore, 
a fully consistent treatment of the perturbations for this model requires to take into account 
perturbations on all these variables at the same time. This is done automatically by the 
third order differential equation (3.34). In sect. |^ we will use this equation for the study of 



the structure formation properties within the running FLRW model. 

4. A FLRW-like model in which only A varies: The A t CDM model 

Let us now consider the possibility of interacting A-cosmology, that is, a cosmology in 
which the vacuum exchanges energy with matter at a fixed value of G = Go. The global 
dynamics of such models has been investigated extensively in the literature on purely 
phenomenological grounds, in fact, much before the discovery of the present accelerating 



stage [21] - see [22| for a review of the old models. However, a variant of these models 



which can be well motivated within the framework of QFT in curved space-time did not 
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appear until later and in subsequent works 2£ ] . Here we briefly present the main 
points of the theoretical analysis - for an updated and thorough discussion, see Basilakos 
et al. I^TJ. Following our general QFT framework of section ||, we expect a functional 
dependence with H, i.e. pa = Pa(H), in which the first non-trivial power with respect to 
the Hubble rate is H 2 , meaning that the basic dependence with H is again of the form (|3.3|), 
or equivalently (3.6). Since G is constant now, the exchange of energy between vacuum 
and matter is governed by equation Q2.9| ). Combining the latter with Friedmann's fl2.4j ) for 
a flat universe and the acceleration law (|2.5| ), one obtains 



H + -H 2 = 4vr G (p A - uj mPr , 



(4.1) 



where oj m = p m / ' p m . In the particular case uj m = (matter dominated epoch), we have f31] 



H+-H 2 
2 



A 



Am G pa . 



(4.2) 



The traditional cosmology (i.e. the concordance ACDM model) can be obtained as a 
particular case of this equation by direct integration of it for p\ =const., but the same 
equation (^) is also valid for p\ = p\(H). Therefore, if we insert ( |3.6| ) in (^) we can 
obtain the corresponding Hubble expansion rate for the A^CDM model. It turns out that, 
in this case, an explicit integration with respect to the cosmic time is feasible (a feature 
which was impossible for the A t G(CDM model of section ||). The final result is |3l]] 



H(t) = H 



coth 



5h„ 

2 



(09- u)(l- u)t 



(4.3) 



From here the corresponding scale factor a{t) can also be found immediately: 
a(t) -- 



n 







3(1-") 



sinh 3 ^-") 



-H y/(SII-U)(l-U)t 



(4.4) 



Eliminating the time variable between equations ( |4.q ) and (|4.4| ) we arrive at an expression 
of the normalized Hubble flow in terms of the scale factor, or equivalently the redshift 
z = (1 — a)/a: 



E{z) 



H(z) 
Hn 



+ 



r 1 + z) 3(l-u) 



nl/2 



(4.5) 



1-v 1 ■ 

Let us finally report on the behavior of the matter and vacuum energy densities in this 
model as a function of the redshift. This is an important point as it shows one of the basic 
differences between this model and the one considered in the previous section. At this 
point we will come back to the general matter EoS parameter uj m = p m /Pm because we 
wish to see the change in the evolution of both cold DM and relativistic particles. Starting 
from the covariant conservation law fl2.9|) , and trading the time derivatives for derivatives 
with respect to the redshift (through the relation d/dt = —H(l + z)d/dz), we obtain 



dfm 

dz 



3(1 + UJmjPr, 
l + z 



dpA 

dz 



(4.6) 
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Using this equation in combination with (|3.q ) and Friedmann's (|2,4j ), we arrive at a simple 
differential equation for the matter density, 

dPm 3(1 +uj m ) (1 - v) 

—. — Pm = , 4.7) 

az 1 + z 

whose trivial integration yields 4 

p m (z) = p° m (l + zf( 1+ ^ 1 -»K (4.8) 

Here p^ m is the matter density at the present time (z = 0). The vacuum energy density 
then follows from integrating once more (|4.6|) using the result Q4.8]): 



p A ( Z ) =pl+ ^fL [(1 + z )3(l+u%0(l-,) _ jj _ (4 9) 

Substituting ( [Q| ) in the last two equations, with z = (1 — a)/a, one may obtain the 
explicit time evolution of the matter and vacuum energy densities, if desired. Recall that 



w m = (u> m = 1/3) for non-relativistic (relativistic) matter. From equation (|4.8|) we 
see that, for the matter epoch, the density does no longer evolve in the standard form 
Pm(z) = Pm(^ + z ) 3 but as p m (z) = p„(l + z) 3 ^ u \ The ^-correction in the power is 
caused by the exchange of energy between matter and the vacuum. This is also reflected in 
the corresponding non-constant behavior of pa( z ) m ( |4.9| ). Similarly, the energy density of 
radiation takes on the corrected form pr(z) = p R (l+z) 4 ( 1 ~ u \ instead of pr(z) = /^(l+z) 4 . 
Such correction to the radiation density can be used to put a moderate limit \u\ < 0.1 to 
the basic parameter v in this model from primordial nucleosynthesis p8[ , although that 
limit is actually superseded by better bounds derived from other precision observables 
which will be discussed in the next section. Obviously, for v = we recover both the 
standard expressions for the evolution of the matter density in the ACDM cosmology, and 
the constancy of the vacuum energy: p\ = p9 . As a consistency check, we note that if we 



substitute (4.9) into (3^6), we recover the normalized Hubble flow (4.5) for this model in 
the matter dominated epoch. 

The following observations are in order. Owing to the coupling between the time- 
dependent vacuum and the matter component in the A^CDM model, there is either a 
particle production process or an increase in the mass of the dark matter particles [^T[ [52| . 
The microscopic details of this coupling are not provided by studying the aforementioned 
global cosmological properties of the model. Therefore, one has to preclude the possibility 
that there can be a significant decay of vacuum energy into matter, especially into ordinary 
cold and relativistic matter (e.g. baryons and photons). This is usually not the case 
because on theoretical grounds the parameter v in ( |3.7| ) is naturally predicted to be small 
(M ^ 10 -2 ) (as it plays the role of a small /3- function) , and moreover in practice the fits 
to the cosmological data also confirm this fact and enforce it to be even smaller. 

If we compare the situation with the AjGtCDM running model of the previous section, 
we can see that both models share the running vacuum law ( |3.6[ ) - and hence the depen- 
dence on the ^-parameter - but in the AjGtCDM case matter is covariantly conserved, and 



4 This equation was first studied in |28j], and later was also considered in |5l| 
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as a result there is no danger that the vacuum decays in excess into ordinary matter. In 
contrast, the feedback associated to the time evolution of the vacuum in the AtGjCDM 
model produces in this case a slow (logarithmic) evolution of the gravitational constant 
with the universe's expansion, see ( |3.11| ). We shall analyze the best fit values of v for 
the two running models A^CDM and A^f^CDM when we compare them with the current 
precision observations in the next sections. 

Finally, as the data on structure formation will be involved in restricting the value of 
let us also quote the corresponding differential equation for the matter density contrast 
D = Spm/Pm in the current A t CDM model. Denoting differentiation with respect to the 



cosmic time by a dot, we have [31] 



D + (2H + Q)D- 4irGp m -2HQ-Q 



D = 0, 



(4.10) 



where the time evolving vacuum energy pa = p\(t) affects the growth factor via the function 



Q(t) = -pa/pt, 



(4.11) 



For pa =const. the equation ( |4.10| ) reduces to the standard form of the matter perturba- 
tions in the ACDM model, i.e. D + 2HD — 4irGp m = [pUfl. Recall that this equation is 
formally valid both for p\ = and for pa 7^ provided H in it is changed accordingly. The 



numerical solution of the more general equation ( 4.10 ) for the case of the A t CDM model 
has been studied in [31], and the corresponding results will be used in the section ^ in 
order to confront the running cosmological models with the observations. Let us point out 
that, for the A^CDM model, the perturbations in pa are subject to some ambiguities ex- 
plained in Ref. [54], and therefore for this model we shall present the matter perturbations 
in the effective framework defined in equation (4.10) , in which it is perfectly consistent to 
have matter perturbations for any non-perturbed vacuum function pa = PA(t)- This is in 
contrast to the A^G^CDM model, where perturbations in all variables are necessary to get 
a consistent picture, as explained in the previous section. 



5. Equation of state analysis of the running vacuum models 

A very important part in the study of the Hubble expansion behavior of a given cosmological 
model is the analysis of the equation of state (EoS) of the dark energy component. In 
contrast to cold and relativistic matter, the effective EoS of the DE may appear in general 
as a non-trivial function of time u; e ff(t), or of the cosmological redshift, uj e s(z). This 
feature may also be expected for general CC models of the form (|2.1|), for which the EoS 



is pa(z) = —pa(z). One can easily understand the reason as follows - see [12] for details. 
The usual way to parametrize the DE is to consider such entity as a self-conserved "fluid" 
at fixed Newton's coupling G. This representation may be called the "DE picture" of the 



original CC model [12]. Since the CC models under consideration do not satisfy at least 
one of these two conditions (either because the vacuum energy interacts with matter or 
because it evolves with time along with G), the effective EoS of the system in the DE 
picture will be time/redshift dependent. 
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Let the self-conserved fluid of the DE picture be characterized by density and pressure 
(p ms ,PD), therefore satisfying p ms + a Hp p ms = and p D + 3 F D (1 + uj eS )p D = 0, where 
u e s{z) = pd(z)/ Pd(z) is the non-trivial EoS function of the redshift that we are looking for. 
The subscript D in these variables serves us to distinguish the DE picture from the original 
CC picture -that is, the original description of the model in terms of variable cosmological 
parameters, which we will denote by H\. Of course p ms ^ p m and po 7^ PA in general. 
The Hubble rate in the DE picture takes on the form: 



Hi 



8irG 



( Pms ( z ) + PD {z)) = Hi n° m (i + z) a + n D (z) 



where 



CId(z) = Q° D exp < 3 



dz 

o 1 + *' 



(5.1) 



(5.2) 



Pms/Pc- Notice that the parameter need not coincide in 



and we have defined (l^ — "° 
general with f2^, as they are determined from two different parameterizations of the data. 
Similarly, we have introduced Qd(z) = pd{z)/p° c - and its current value Sl° D = Qd(0). 
Notice also that Go = 1/Mj, is strictly constant, whereas G in (|3.16| ) is in general a 
variable function of the redshift: G = G(z). The effective EoS in the DE picture now 
follows from 

1 1 + z dpD 



-1 + 



(5.3) 



3 pn dz 

where it is to be understood that po must be computed from the matching condition 
between the DE picture and the original CC one, i.e. by requiring that the expansion 
histories of the universe are numerically equal in both pictures: H-q{z) = H\{z) . Using 
the above formulas, let us compute the effective EoS in each case for the matter dominated 
epoch. For the running FLRW model (A^GfCDM) of section pi we find 



(g(z)n° m - n° m ) (i + z) 



(5.4) 



where 0,\(z) is given by ( p. 25 ) and g(z) by ( p. 28 ), and use has been made of the Bianchi 
identity ( |3 . 1 9| ) . Similarly, for the running AjCDM model of section ^j, we obtain 



^ S {z) = -1 + (1 



n° m (i + zf^) 



1 + z) 



n° m [(i + ^a-") -i] 



1 



:i+z) 



1] 



(5.5) 



In these formulas, we assume that the cosmological mass parameters in the two pictures 
are in general different, i.e. AQ.^ = 0,^ — ^ 0. 

It is remarkable to note that in all time-dependent vacuum models of the form (^[l]) , 
a value z = z* always exists for which the effective EoS ui e g(z) crosses the phantom divide 
u e fi(z*) = —1, see [^] for the general proof of this statement. For the models under 
consideration we can easily check that it is so. This is particularly simple if we look at 
equation (|5.4j ). For, as 0,^ and must be very close, and g(z) is a slowly monotonous 
function, the effective EoS crosses the phantom divide at some point z* characterized by the 
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Figure 1: The behavior of the effective EoS function u> e g(z) for the two running vacuum models. 
In each case, corresponds to the central value obtained in our fit (cf. Sect. whereas the value 
of v is in the lower limit of the la range. For fixed v — —0.004, and different values and signs of 
Af2 m , the effective EoS behavior of both running models appears as significantly different from the 
ACDM prediction (ua = — 1). Indeed, the effective EoS of the running A models is seen to present 
quintessence- or phantom-like behavior, depending on the values of the parameters. 



condition = O^. Whether the crossing point z* is in the past or in the future, it 

will also depend on the sign of AS7^ and the value of v. Clearly, for AO^ = the crossing 
is exactly at z* = since g(0) = 1. Similarly, it is easy to see that, for z — > 0, equation 
J5.5Q renders u>^ s — > — 1 for AS7^ = 0, whereas if Af2^ > or Af2^ < the behavior is 
w eff ~ ~~ 1 (quintessence- like) or uj^ s < — 1 (phantom-like) respectively. Furthermore, for 
both models one can easily show that, if AJ7^ = 0, the quintessence/phantom-like behavior 
at points z near our time is completely controlled by the sign of v. Indeed, expanding the 
above expressions in the aforementioned conditions and for small v, we meet: 

co^ G (z) = -1 - 3^ (1 + zf In (0^(1 + zf + fli) + 0(v 2 ) , (5.6) 
no 

u, e A ff (z) = -1 - 3u-^ (1 + zf ln(l + z) + 0{v 2 ) . (5.7) 

These two equations are very similar, but not quite. Nevertheless we find once more that 
both running models share the property w e g = —1 at z = 0. Furthermore, for both 
running models we also learn that, if Af2^ = 0, the cases v < or v > amount to 
quintessence-like or phantom-like behaviors respectively. For fixed v, the quintessence-like 
or phantom-like behavior can also be modulated by the sign of AQ m . Some numerical 
examples are considered in Fig. |], where the crossing is in the future, except for the cases 
AO^ t = in which the crossing is at z = (i.e. precisely at our time). 
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The previous discussion clearly illustrates, with the help of the two non-trivial running 
models AjGjCDM and A 4 CDM discussed in the previous sections, that a generic cosmolog- 
ical model ( |2.lD with time- varying vacuum energy generally leads, when described in the 
DE picture (i.e. as if it consisted of a cosmological self-conserved DE fluid at constant G), 
to a crossing of the phantom divide. Thus it may provide a natural explanation for the 
observational data, which still admit a tilt in the phantom domain ||, or it may simply ex- 
plain the apparent quintessence-like behavior, despite there is no fundamental quintessence 
or phantom field in our framework. In fact, remember that the original CC picture for both 
models is just a time-varying vacuum model with or without interaction with matter. 

In the next section, we confront these running models with the observational data on 
SNIa, CMB and BAO. As we will see, for values of v sufficiently large within the la range 
of our fits, both models predict a cluster redshift distribution that differs noticeably from 
the ACDM expectations. For such values of u, the EoS function of our models should also 
be easily distinguishable from the ACDM value, ui\ = —1, even in the simplest = 

case. We refer once more to Fig. [l], where we have plotted the EoS functions ( |5.4j ) and 
( |5.5| ). For large redshifts [z ~ 1.5 — 2) it is even conceivable that the quality of the data 
provided by forthcoming SNIa experiments could be enough so as to discriminate between 
both running cosmological models. A combined analysis of clusters and the dark energy 
EoS could be a most effective strategy to put our models to the test. 



6. Confronting the running vacuum models with the latest observations 

In this section, we briefly present the basic observational samples and data statistical 
analysis that will be adopted to constrain the "running" models presented in the previous 
sections. First of all, we use the Constitution set of 397 type la supernovae of Hicken et 
al. |||. In order to avoid possible problems related with the local bulk flow, we use a 
subsample of 366 SNIa, excluding those with z < 0.02. The corresponding XsNia function, 
to be minimized, is: 



366 

XSNIa(P) = Yl 



M th KP)-M° bS (V n2 



(6.1) 

= i L *i 1 

where ai = (1 + Zj) _1 is the observed scale factor of the Universe for each data, zi is the 
observed redshift, fi is the distance modulus [i = m — M = 51og<i£ + 25 and oIl(cl, p) is the 
luminosity distance: 

dL ^ p) = -aI a (6 - 2) 

with c the speed of light and p a vector containing the cosmological parameters (i.e., f^, v) 
that we wish to fit for. The previous formula applies only for spatially flat universes, which 
we are assuming throughout. 

On the other hand, we consider the baryonic acoustic oscillations (BAOs), which are 
produced by pressure (acoustic) waves in the photon-baryon plasma of the early universe, 
as a result of the existence of dark matter overdensities. Evidence of this excess has been 
found in the clustering properties of the SDSS galaxies (see [j55f| , f56|| ) and it provides a 
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"standard ruler" we can employ to constrain dark energy models. In this work we use the 
newest results of Percival et al. @, r s (z d ) / D w (z±) = 0.1390 ± 0.0037 (see also §7|]). Note 
that r s is the comoving sound horizon size at the baryon drag epoch 5 zj, Dy(z) is the 
effective distance measure [55| and z+ = 0.275. Of course, the quantities (r s ,Dy) can be 
defined analytically. In particular, r s is given by: 

" ad da 



r s(zd) = —e- 



c 

V3 



a 2 H{a)Jl + (3ft0/4f}0) a 



(6.3) 



where = (1 + Zd) 1 and VL®h? ~ 0.02263. In this context, the effective distance is (see 
H): 

1/3 



Dy(z) 



(l + zYDi(z) 



cz 
Wz) 



(6.4) 



where Da{z) = (1 + z) 2 di(z,p) is the angular diameter distance. Therefore, the corre- 
sponding Xbao function is simply written as: 

1 2 



Xbao(p) 



0.1390 



0.0037 2 



(6.5) 



Finally, a very accurate and deep geometrical probe of dark energy is the angular scale 
of the sound horizon at the last scattering surface, as encoded in the location lf T of the first 
peak of the Cosmic Microwave Background (CMB) temperature perturbation spectrum. 
This probe is described by the CMB shift parameter ]59| , |6Q ], defined as: 

„ — f 1 da , 

R=V^ 6.6 
Ja ls a 2 E(a) 

The measured shift parameter according to the WMAP 7- years data || is R = 1.726±0.018 
at z ls = 1091.36 [or a ls = (1 + z^)- 1 ~ 9.154 x 10" 4 ]. In this case, the X cm b function is 
given by: 

2 [R(p) - 1.726] 2 

Xcmb(P) = • (6-7) 

Note that the measured CMB shift parameter is somewhat model dependent, but mainly 
for models that include massive neutrinos or those with a strongly varying equation of state 
parameter (which is not our case). For a detailed discussion of the shift parameter as a 
cosmological probe, see |pj| . 

In order to place tighter constraints on the corresponding parameter space of our 
model, the probes described above must be combined through a joint likelihood analysis 6 , 
given by the product of the individual likelihoods according to: 

£tot(p) = ^SNIa X ^BAO X £cmb , (6.8) 



5 Zd is given by the fitting formula of |5q ]. 

6 Likelihoods are normalized to their maximum values. In the present analysis we always report la 
uncertainties on the fitted parameters. Note also that the total number of data points used here is N to t = 
368, while the associated degrees of freedom is: dof= Ntot — nut, where nat is the model-dependent number 
of fitted parameters. 
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Figure 2: Upper Panel: Likelihood contours (for — 21n£/£ max equal to 2.30, 6.16 and 11.81, 
corresponding to la, 2a and 3a confidence levels) in the (f2° n ,^) plane for the A t GtCDM model. 
The left panel shows the contours based on the SNIa data (thick solid black lines), BAOs (dotted- red 
lines) and CMB shift parameter (dashed-green lines) . In the right panel we show the corresponding 
contours based on the joint statistical analysis (SNIa+BAO+CMB data). Lower Panel: Similar 
likelihood contours for the A t CDM model. 



which translates in an addition for the joint x 2 function: 

XL(P) = XSNIa + XBAO + Xcmb ■ ( 6 - 9 ) 

In our x 2 minimization procedure, we use the following range and steps of the fitted 
parameters: G [0.1, 1] in steps of 0.001 and v G [-0.02,0.02] in steps of 10" 4 . 

In Figure |2| we present the results of our analysis for the A(G(CDM (upper panel) and 
the A^CDM (lower panel) models. The left plots in that figure show the la, 2a and 3a 
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confidence levels in the (O^, v) plane, with the SNIa-based results indicated by thick solid 
lines, the BAO results by dotted - red lines and those based on the CMB shift parameter by 
dashed -green lines. Using the SNIa data alone it is evident that although the parameter 
is tightly constrained (~ 0.284), the v parameter remains completely unconstrained. As 
can be seen in the right plots of Figure g, the ab ove degeneracy is broken when using the 
joint likelihood analysis, involving all the cosmological data. For the A t G^CDM we find 
that the overall likelihood function peaks at = 0.283+o;o}g and v = -0.001 ± 0.0032 
with Xtoti^m^) — 439.36 for 366 degrees of freedom. 

On the other hand, for the AjCDM model the best fit parameters are: = 0.284^q'q}4 
and v = -0.0017^;^3 with Xtoti^m, v)/dof ~ 439.34/366. Finally, for the usual A 
cosmology (v = and G = Gq) we find that the overall likelihood function peaks at 
n° m = 0.284 ± 0.012 with x1ot(^m)/ do f - 439.61/367. We remark that although for both 
running models the best fit value for the v parameter is negative 7 , its la range includes 
both positive and negative v values and encompasses quite a wide segment. For instance, 
for the A^GfCDM model the la range is roughly [—0.004,0.002]. It is worth mentioning 
that this range of values is in the same ballpark as the bounds obtained for v from the 
normalization of the spectrum amplitude and from primordial nucleosynthesis |l^] . 

To end this section, let us mention that we have checked that using the earlier BAO 
results of Eisenstein et al. [55| does not change significantly the previously presented 
constraints. 



7. Halo abundances and their evolution in dark energy models 

In an attempt to define observational criteria that will enable us to identify the realistic 
variants of the concordance ACDM cosmology and distinguish among the different A(t) 
models, we derive and compare their theoretically predicted cluster-size halo redshift dis- 
tributions. 

The formalism to compute the fraction of matter in the universe that has formed 
bounded structures and its redshift distribution was developed in 1974 by Press and 
Schechter |32| (hereafter PSc). In their approach, the primordial density fluctuation for 
a given mass M of the dark matter fluid is described by a random Gaussian field. One 
introduces a function, J-(M, z), representing the fraction of the universe that has collapsed 
by the redshift z in halos above some mass M. With this function one may estimate the 
(comoving) number density of halos, n(M, z), with masses within the range (M, M + 6M): 

n(M,z)dM = ^^ldM. (7.1) 
This expression can be rewritten as follows: 



7 Notice that in Ref.lpq] the best fit value for the A t CDM model was v ~ 0.002 > 0, because v was 
actually sampled only within the interval v G [0, 0.3]. The value of fi^i however, remains very similar. 
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where /psc("") = y/^/^i^c/^) exp(— 5 2 /2a 2 ), 5 C is the linearly extrapolated density thresh- 
old above which structures collapse ||| (see Appendix [A| for more details), while a 2 (M, z) 
is the mass variance of the smoothed linear density field, which depends on the redshift z 
at which the halos are identified. It is given in Fourier space by: 

a\M,z) = I k 2 P{k)W 2 {kR)dk , (7.3) 

2vr 2 J 

where D(z) is the growth factor of perturbations 8 , P(k) is the power-spectrum of the 
linear density field, and W(kR) = 3(sinkR — kRcoskR) / (kR) 3 is the top-hat smoothing 
function, which contains on average a mass M within a radius R = (3M/47T/9) 1 / 3 and 
p = 2.78 x lO 11 fi° l /t 2 M Mpc -3 . We use the CDM power spectrum: 

P(k) = P k n T 2 (n° m ,k), (7.4) 

with T(n° m ,k) the BBKS transfer function [§§]: 

k) = ln(1 2 + 3 ^ 34g) [l + 3.89, + (16.1g) 2 + (5.46,) 3 + (6.71g) 4 ] " V * , 

q = q(k) = A . (7 .5) 

Here T is the shape parameter, according to [65], provided by the 7-year WMAP results @, 
in combination with the model fitted f2 m values of section 6. Note that in this approach 
all the mass is locked inside halos, according to the normalization constraint: 

/+oo 
/psc^dlmr- 1 = 1 . (7.6) 
-oo 

It is traditional to parametrize the mass variance in terms of the rms mass fluctuation 
amplitude on scales of Rs = 8 h~ l Mpc at redshift z = [ag = as(0)}, so that we have: 

2 (M , J °°k^ 2 T 2 (nl,k)W 2 (kR)dk 

where 

D(z) . . 

mz) = a sjj^ ■ ( 7 - 8 ) 

Although the previously described Press-Schechter formalism was shown to provide a 
good first approximation to the halo mass function obtained by numerical simulations, it 
was later found to over-predict /under-predict the number of low/high mass halos at the 



present epoch [|66| , 67 1. More recently, a large number of works have provided better fitting 



The equations to compute the growth factor have been given at the end of Sects. ^ and namely 



Eqs. (3.34) and (4.10). Since the pure matter universe (Einstein de-Sitter) has the solution D^dS = a , we 
normalize our vacuum models such as to get D ~ a at large redshifts (e.g. z — 500), where the matter 
component dominates the cosmic fluid and our models resemble a pure CDM model. 
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functions for f(cr), some of them based on a phenomenological approach. In the present 
treatment, we adopt the one proposed by Reed et al. 



f{cr,n eS ) = A\ — 



a 2 V 

1 + ( ^| J + O.6G1 + 0.4G 2 



5c 

— exp 
a 



0.03 



ebb 2 . 



la 2 (n eff + 3) 2 \a 



0.6' 



(7.9) 

where A = 0.3222, p = 0.3, b = 0.707, c = 1.08, while Gi, G 2 and n eff , the slope of the 
non-linear power-spectrum at the halo scale, are given by: 



G\ = exp 



(In a' 1 - 0.4) 5 
2(06)2 



, G 2 = exp 



(In a" 1 -0.75) 2 
2(0^)2 



diner" 1 
"eff = 6— — — 3. 



dlnAf 



(7.10) 



7.1 Collapse threshold and mass variance of the running vacuum models 

In order to compare the mass function predictions for the different vacuum models, it is 
imperative to use for each model the appropriate value of 5 C and a$. Indeed, the way in 
which dark energy (time-varying vacuum energy, in our case) affects the formation of gra- 
vitationally bound systems (clusters of galaxies) is a crucial question that has received a lot 
of attention in the literature. Three distinct scenarios (with increasing level of difficulty in 
the theoretical treatment) can be considered, to wit: (i) dark energy remains homogeneous 
and only matter virializes; (ii) dark energy may cluster but only matter virializes; and (iii) 
dark energy may also cluster and the whole system (matter and dark energy) virializes. 
The first scenario has been widely used in the literature, probably because it is the simplest 
one [35, 



37, 38 



|3J|. In the present paper, we will consider just this canonical scenario 
when dealing with the AfCDM model. For the full running FLRW model (or A^G^CDM 
model), however, which is the main object of the present study, we will go a bit deeper 
and we shall treat the gravitational collapse within the framework of the second scenario 
(in which SG and 5\ are both non- vanishing and only matter virializes). Then, in the 
Appendix ^ we shall compare these results for the A^G^CDM model within scenario (ii) 
to those obtained within the canonical scenario (i), for which 5G = 5a = 9 . 

It is well known that for the usual A cosmology 5 C ~ 1.675, while Weinberg & 



Kamionkowski [34] provide an accurate fitting formula to estimate 5 C for DE models with 
constant EoS parameter. It should furthermore be noted that such a formula is, in princi- 
ple, only applicable to "conventional" DE models, i.e. models with self-conserved matter 
and DE densities and constant G. Although these conditions are not satisfied by our run- 
ning vacuum models, one can show [38, [!(]] that the 5 C values for a large family of dark 
energy models with a time-varying EoS parameter can be well approximated using the 
previously discussed fitting formula, as long as the EoS parameter is not very different 



9 A complete study of all the possible clustering scenarios in the context of our vacuum models is beyond 
the scope of the present paper. Scenario (iii) is of course the most difficult one, and although a first attempt 
to include the effects of clustered and globally virialized vacuum energy and matter for the AtCDM model 
has been considered in pH , the modeling of this situation requires the introduction of new assumptions 
and more free parameters. In order to avoid a too cumbersome treatment here, the corresponding study 
for the AtGtCDM model will be presented elsewhere. 
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from -1 near the present epoch. This is indeed the case for our models, as shown in Fig. [T]. 
In spite of this prognosis, we will nevertheless carefully compute the value of 5 C within our 
vacuum models, following the prescriptions given in [|3S( | (cf. Sect. 2.1 in that reference) 
and [33]. Starting from a Newtonian formalism, we derive a non- linear second-order differ- 



ential equation for the evolution of the matter perturbations in our vacuum models from 
which we compute S c (z), and in particular S c = 5 C (0) (cf. Appendix |A|for the calculational 
details and methodology). The resulting values are listed in Table 1. As an example, in 
the case of the A t G t CBM model, for £l° m = 0.283 and v = -0.001 we find 5 C ~ 1.677, 
pretty close to the ACDM value (~ 1.675). On the other hand, for the A^CDM model 
(0°, = 0.284, v = -0.0017) we find a bit larger value, 5 C ~ 1.685. In all cases, however, 
the differences with respect to the ACDM value are at the few per mil level only. This is in 
agreement with the expectations for our models, and it is exactly the same that happens 



in other models [35, 40]. 



The relevant as value for the different "running" vacuum models can be estimated by 
scaling the present time ACDM value 10 (o"s,a) using equation (7^3). In general, for any DE 
model, equation ( [7.3] ) takes the form: 

1/2 



&8 = <?8,A- 



P J °° k n+2 T 2 (n° m , k)W 2 (kR 8 )dk 



Pa,o / °° kn+2T2 ^°m,^ k)W 2 (kR 8 )dk 



(7.11) 



•£>a(0) 

where P /P Afl = « A /^) 2 . 

Therefore, using the observationally determined value for cg^, we can easily derive 
the corresponding <jg value for any of the time-varying vacuum models. Indeed, the joint 
WMAP7+BAO+i?o analysis of Komatsu et al. || resulted in a value of cr 8)A = 0.811. This 
value is in quite good agreement with that provided by a variety of different methods; for 
example, an analysis based on cluster abundances gave the degenerate combination: a 8) \ = 
(0.83 ± 0.03) (n^/0.25) -0 - 41 ||. The weak-lensing analysis of Fu et al. @ provided 
0"8,A = (0.837 ± 0.084) (r^/0.25) -0 ' 53 , while two recent studies based on a joint analysis 
of the large-scale clustering of red SDSS galaxies, CMB, SNIa and BAO data resulted in: 
cr 8) A ^ 0.8 ± 0.02 (for 0°„ ~ 0.26) fn], [T2|. The only method that provides discrepant 



results is that based on galaxy or cluster peculiar velocities [|73|, |74j] . 

Finally, inserting the Komatsu et al. 0] ag,A value in equation ([7.1 1|) we can estimate 
the corresponding as values for our models. These values will be used in the mass function 
analysis and are listed in Table 1. 

7.2 Halo mass function & number counts of the running vacuum models 

We now move to our results. Given the halo mass function from Eq. ( [7.2| ) we can derive 
an observable quantity which is the redshift distribution of clusters, M(z), within some 
determined mass range, say Mi < M < Mi- This can be estimated by integrating the 
expected differential halo mass function, n(M,z), with respect to mass, according to: 



dV 



Ah 



M{z) = — / n(M,z)dM, (7.12) 
dz J Mi 

10 In the following discussion, the quantities referred to the ACDM model are distinguished by the subscript 
'A' (erg, a; -Da; Hi a) whereas the corresponding quantities in the running models carry no subscript. 
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where dV/dz is the comoving volume element, which in a flat universe takes the form: 

dV n, . driz) 

4vrr 2 (z)-^, (7.13) 
dz dz 

with r(z) denoting the comoving radial distance out to redshift z: 

. . c f z dz' , w i A . 

^> = wj w-y (7 ' 14) 

In the upper panel of Fig. ||, we show the theoretically expected redshift distribution, N(z), 
for cluster-size halos, ie., M 1 = 10 134 h~ l M Q and M 2 = 10 16 h~ l M & for the A^CDM, 
A^CDM and ACDM models, using the best-fit values for the and v parameters provided 
in section 6 (also indicated in Table 1). In the lower panel we show the relative differences 
of the two running models with respect to the concordance ACDM model. The different 
models are characterized by the symbols and line types presented in Table 1. 

It is evident that, for the central fit values of the v parameter, the A^CDM is the 
only model which shows significant differences with respect to the concordance ACDM 
model, producing a larger number of cluster-size halos at any given redshift. As for the 
AtG(CDM model, although it also produces larger numbers of halos than the ACDM model, 
the differences are minimal. For instance, for 2 < z < 3 the A^CDM model presents 
8N /N ~ 0.47, whereas in the A^GjCDM case the corresponding value amounts to only 
~ 0.06. 

However, in order to assess whether it is possible to observationally discriminate be- 
tween the different vacuum models, we need to concentrate on a more significant model 
comparison by deriving the expected halo redshift distributions in the context of realistic 
future cluster surveys. Two of these realistic future surveys are: 

(a) The eRDSITA satellite X-ray survey, with a flux limit of f\ im = 3.3 x 10~ 14 ergs s~ 4 
cm -2 , at the energy band 0.5-5 keV and covering ~ 20000 deg 2 of the sky. 

(b) The South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) survey, with a limiting flux 
density at vq = 150 GHz of f uo ,\\ni = 5 mJy and a sky coverage of ~ 4000 deg 2 . 

To realize the predictions of the first survey we use the relation between halo mass and 
bolometric X-ray luminosity, as a function of redshift, provided in Jf5[: 



L(M,z) = 3.087 x 10 44 



ME(z) 



1.554 



hT 2 ergs -1 . (7.15) 



10 15 h- l M Q _ 

The limiting halo mass that can be observed at redshift z is then found by inserting in the 
above equation the limiting luminosity, given by: L = ^irdj^fi^Cb, with d^ the luminosity 
distance corresponding to the redshift z and q, the band correction, necessary to convert 
the bolometric luminosity of Eq. (7.15) to the 0.5-5 keV band of eRDSITA. We estimate this 
correction by assuming a Raymond-Smith |7^] plasma model with a metallicity of 0AZ & , 
a typical cluster temperature of ~ 4 keV and a Galactic absorption column density of 
tih = 10 21 cm -2 . 

The predictions of the second survey can be realized using again the relation between 
limiting flux and halo mass from [75|: 

2.592 x 10 8 mJy ( M V' 876 f2/3 ^ f - ^ 

^ = — dW) — V^m~J {z) { ] 
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Figure 3: The theoretical redshift distribution of MJ> 10 13A h~ 1 M Q clusters (upper panel) for 
the different vacuum models and the corresponding fractional difference (lower panel) between the 
running models and the reference ACDM model. The running models are plotted for the central fit 
values of the v parameter, v = —0.0017 (red dashed line) for the A 4 CDM and v = —0.001 (black 
crosses) for the AjGtCDM. The full list of input parameters, corresponding to the central fit values, 
as well as the symbols or line- types characterizing the different models, are shown in Table 1. Error 
bars are la Poisson uncertainties. 



where d_&(z) = di/{\ + z) 2 is the angular diameter distance out to redshift z. 

In Fig^ (upper panels) we present the expected redshift distributions above a limiting 
halo mass, which is M\ = Mi; m it = max[lO 13 ' 4 /i _1 M0, Mf], with Mf corresponding to the 
mass related to the flux- limit at the different redshifts, estimated by solving Eq. ( 7.15 ) and 
Eq. ( [7.16 ) for M. In the lower panels we present the fractional difference between the two 
time-varying vacuum models (A^CDM and A^G^CDM) and the ACDM, similarly to Fig.|||, 
but now for the realistic case of the previously mentioned future cluster surveys. The error 




Figure 4: The expected cluster redshift distribution, over the whole sky, of the two running 
vacuum models for the case of the two future cluster surveys eROSITA and SPT (upper panels), and 
the corresponding fractional difference with respect to the reference ACDM model (lower panels). 
The input parameters, corresponding to the central fit values, as well as the symbols or line-types 
characterizing the different models are shown in Table 1. 



bars shown correspond to 2a Poisson uncertainties, which however do not include cosmic 
variance and possible observational systematic uncertainties, that would further increase 
the relevant variance. 

It is evident that the imposed flux limits, together with the scarcity of high-mass halos 
at large redshifts, induces an abrupt decline of M{z) with z, especially in the case of the 
eROSITA X-ray survey (note the shallower redshifts depicted in Fig.|| with respect to Fig.||). 
Furthermore, in the case of the A^CDM model the uncertainties indicate that we should 
expect significant and large differences to be measured in both experiments for any redshift 
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Figure 5: Extension of the analysis presented in Fig. 4, but for the A t GtCDM model only. Once 
more the theoretical redshift distribution of M<; 10 13 Ji^Mq clusters in the two future surveys 
is shown in the upper panel, whereas the corresponding fractional difference with respect to the 
ACDM model is shown in the lower panel. In this case, we use the limiting values of v comprised 
within the la range of our fit, namely v = —0.004 (black squares) and v = +0.002 (green triangles). 
See Table 1. 



<; 0.3-0.4. On the other hand, the deviations of the AjG^CDM model with respect to the 
concordance model, although significant, are very small (for the central fit values) and thus 
more difficult to detect. 

Let us stress that the previous analysis of the cluster-size halo redshift distribution 
has focused on the central values of our fits to the SNIa+BAO+CMB data (cf. section 
|6|). However, there is no compelling reason to prefer the central values to any other value 
within the la contour. This observation is particularly important when the dispersion 
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Model symbol 




(57V/7V A )cROSITA 

z < 0.2 0.5 < z < 1 


(5JV/M a )spt 
z < 0.2 1 < z < 2 


ACDM black line 
Af G^CDM black crosses 
AfCDM red dashed line 


0.284 0.811 1.675 
0.283 -0.001 0.818 1.677 
0.284 -0.0017 0.840 1.685 


0.00 0.00 
0.01 0.04 
0.05 0.25 


0.00 0.00 
0.01 0.05 
0.06 0.36 


A^G^CDM black squares 
AfGjCDM green triangles 


0.283 -0.004 0.853 1.683 
0.283 0.002 0.786 1.671 


0.06 0.32 
-0.05 -0.18 


0.07 0.45 
-0.05 -0.23 



Table 1: Numerical results. The 1 st column indicates the dark energy model. The 2 nd column 
shows the symbols or lines used to represent the different models in Figs. [| |] and |[ The 3 rd , 
4 th and 5 th columns list the central values of the input fitted cosmological parameters (first three 
rows, one for each model being considered) and their limiting la values (last two rows, only for the 
A t G f CDM model) obtained by using the SNIa+BAO+CMB data (cf. section g). The procedure 
to compute the critical overdensity 6 C for each model (6 th column) is explained in Appendix |A|. 
The remaining columns present the overall, within the indicated redshift bins, fractional relative 
difference [SAF/AF A = (A/ - — A/a) /A/a] between the cluster abundances predicted by the time varying 
vacuum models and by the ACDM cosmology, for the two future cluster surveys discussed in the 
text. No errors are shown since the 2a Poisson uncertainties are less than 10 -2 . 



turns out to be considerably large in that range. As we saw in section 6, the la range 
of the fitted v parameter for both running models is indeed quite large and includes both 
positive and negative values. In particular, for the A^GfCDM model we have obtained v G 
[—0.004,0.002]. This model deserves special attention in that, in contrast to the A^CDM, 
the time evolution of the vacuum energy density is compatible with matter conservation (a 
feature that shares with the standard ACDM model). It is therefore worthwhile to perform 
a more exhaustive investigation of the cluster halo redshift distribution for this model 
within the limits of the la range, so as to check if the predicted deviations with respect to 
the ACDM can be significant. We display the results of this extended analysis in Fig. [5], 
where we show, following a pattern entirely similar to that in Fig. ||, the expectations 
of the A^GjCDM model for the limiting v values in its la range [ie., v = —0.004 (black 
squares) and v = 0.002 (open triangles)]. It is evident that, in this case, the predicted 
redshift distribution of clusters, Af(z), for the A^GtCDM model does show significantly 
larger, and potentially measurable, differences (~ 50%) with respect to the concordance 
ACDM model, even at moderate redshifts. 

In Table 1, we show a more compact presentation of our results, including the relative 
fractional difference between the "running" vacuum models and the ACDM model, in two 
characteristic redshift bins, and for both future surveys under consideration (eROSITA and 
SPT). The first redshift bin corresponds to the local universe (z<? 0.2) and it is shown in 
order to allow us to judge whether the existing cluster samples could be used to discriminate 
among the different models. It is evident that in the local range z<^ 0.2, although there are 
deviations among the models, their amplitude is very small and it will be rather difficult to 
detect them unambiguously. However, at larger redshifts, the amplitude of the deviations is 
such that it should be relatively straightforward to distinguish among the different models. 
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To conclude, in view of the results of our analysis of the cluster halo redshift distribu- 
tion, Af(z), presented in Figures ^, || and Table 1, we can say that both running models 
show large and significant differences with respect to the ACDM within the la range of 
the fitted values of the fundamental parameter v of these models. We have demonstrated 
that these deviations could be measured for the two future cluster surveys considered in 
our study (eRDSITA and SPT) provided we focus on sufficiently high cosmological redshifts 
z^ 0.3, therefore outside the local domain. In practice, eROSITA will be most efficient in 
the redshift range 0.3 < z < 1.2 whereas SPT in the wider segment 0.3 < z < 2.2, where in 
both cases the upper limit on z is determined by the sharp decrease in the cluster number 
counts and the corresponding larger statistical error. 

Let us finally point out that, had we used the same rms mass fluctuation normalization 
for both time- varying vacuum models, i.e., a$ = as,A (as in pHl), we would have found that 
none of the time- varying vacuum models could be distinguished from the reference ACDM 
model at a significant level. This implies that the value of erg plays also a fundamental 
role in these kind of studies, especially when the models have similar Hubble functions and 
perturbation growing modes, as it is the case in the presently studied time- varying vacuum 
models. Therefore, our including the computation of the value of as for each of the models 
under consideration has been an important feature of our study, and illustrates the general 
need to follow this practice for other cosmological models. 

8. Conclusions 

In this work, we have analyzed the observational status of the "running FLRW model" 
(denoted by AfGtCDM), i.e. the cosmological model characterized by the vacuum energy 
density evolving quadratically with the expansion rate, pa(H) = n + ri2 H 2 (with uq ^ 0). 
The non-vanishing coefficient provides the time evolution of the cosmological term, 
and since it is a dimensional quantity it is conveniently parametrized as n<i = (3z^/87r) Mp, 
where v plays the role of the dimensionless /3-function of the running p\. This is particularly 
clear also from the simultaneous running of the gravitational coupling in this model, which 
evolves logarithmically with the expansion rate according to G~ l {H) ~ Mj>(l + ulnH 2 ). 
These two running laws are intimately correlated such that the Bianchi identity is auto- 
matically satisfied, as demanded by general covariance. As a result, a nice feature of this 
model is that matter is locally and covariantly conserved, similarly to the standard ACDM 
model. It means that, in the A^GtCDM model, there is no decay of vacuum energy into 
matter or vice versa. This is in contradistinction to a previous version of the model, the 
running A^CDM model, which was recently confronted with the latest observations in pl| . 
We have re-analysed this model here, together with the AjG(CDM model, so as to better 
highlight the similarities as well as the important differences between them. In both run- 
ning models, Ph{H) is varying formally in the same way as a function of H, but in the 
A^CDM model G is strictly constant at the price of permitting a continuous exchange of 
energy between vacuum and matter. 

We have confronted A t GjCDM and A t CDM in the light of the latest high-quality 
cosmological data from distant type la supernovae, baryonic acoustic oscillations and the 
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cosmic microwave background anisotropics. This has allowed us to put a limit on the size of 
the fundamental parameter v controlling the running of the vacuum energy in both models. 
This parameter also controls the running of the gravitational coupling in the A t GtCDM 
case, or alternatively the exchange of vacuum energy and matter in the A t CDM model. 
Despite the fact that the two models are qualitatively quite different, for both of them 
we find that if v < the formation of structure is reinforced with respect to the ACDM, 
whereas for v > it is depleted. This is understandable in the sense that, for v < 0, 
the vacuum energy decreases in the past (it may even become negative) and therefore the 
formation of structure is favored. In contrast, for v > the vacuum energy becomes larger 
and positive in the past, preventing the growth of structure. 

Although the central fitted v values of the two running models are different, in both 
cases we meet the preferred sign v < 0, with a magnitude of order \v\ = 0(1O~ 3 ). Fur- 
thermore, the corresponding la ranges are quite consistent among them, especially when 
one takes into account their relative uncertainty, which is of the same order as the central 
value, i.e. \Su\ = O(10~ 3 ). As a result, both the A t CDM and A^CDM models could 
show significant departures with respect to the ACDM power spectrum within one stan- 
dard deviation of the fitted values. This fact translates then into a measurable impact on 
the redshift distribution of cluster-size halos, as we have verified in detail, finding that both 
running models could lead to very important deviations with respect to the concordance 
model. 

Interestingly enough, the redshift distribution of cluster-size halos will be measured by 
two future important surveys, one is based on the X-ray eROSITA satellite and the other 
on the Sunyaev-Zeldovich observations performed by the South Pole Telescope (SPT). 
Our analysis shows that by sampling v within its la range there is a significant maximal 
deviation -positive or negative, depending on the sign of v— in the predicted abundance of 
clusters with respect to the ACDM, which could amount to an anomaly of 50% or higher 
at redshifts ranges (z > 0.3) which should be perfectly accessible to both realistic surveys. 

Finally, let us emphasize that both models display an effective equation of state be- 
havior which can mimic both quintessence and phantom dark energy, without of course 
involving quintessence or phantom scalar fields. It follows that the sole EoS behavior of 
the running cosmologies can be highly distinctive with respect to the concordance ACDM 
model and it can also be used to distinguish between the two running models themselves. 
This strategy should be most efficient upon combining the effective EoS determination with 
the predicted deviations of the clustering redshift distribution with respect to the ACDM. 
The upshot of our study is that the "running cosmologies" could provide an alternative and 
successful version of dynamical dark energy which should be testable in the next generation 
of cosmological experiments. 
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A. The critical overdensity in time- varying vacuum models 

In this appendix we explain in detail the computation of 5 C (the linearly extrapolated 
density threshold above which structures collapse) for the RG models we have considered 11 . 
The quantity 5 C is used for the study of the halo abundances and their evolution in Sect.[7|. 

The A t G t CDM model: 

First we have to derive the second-order differential equation governing the non-linear 
evolution of the matter perturbations in our model. Using the Newtonian formalism for 
the cosmological fluid, we start by writing the continuity, Euler and Poisson equations in 
the matter dominated epoch: 

^ +Vr ,.( Pm v) = 0, (A.l) 

dv 

— + (v-V r )v + V^ = 0, (A.2) 
V 2 $ = 47^^^(1 + 3^) , (A.3) 

i 

where v is the total velocity of the co- moving observer in three-space, $ is the Newtonian 
gravitational potential, r is the physical coordinate, Gn is the Newton's coupling and Y2i 
runs over all the energy components, in our case non-relativistic matter and the running 
cosmological constant (i = m, A). Finally, 0Ji = Pi/pi is the EoS parameter for each 
component (0 for dust, and —1 for the CC, respectively). Let us recall that, within this 
framework, Gjy and p\ depend on time (at the background level). We introduce comoving 
coordinates x = r/a and define the perturbations in the following way: 

Pi (x,t) = pi(t) + 6pi(x,t) = pi(t)(l + 5i(x,t)) , (A.4) 

$(x,t) = $ (x,t) + <p(x,t) , (A.5) 

v(x,t) = a(t)[H(t)x + u(x,t)] , (A. 6) 

G N (x, t) = G(t) + SG(x, t) . (A.7) 



11 We follow the standard methods available in the literature - see e.g. |3q] and j33|], and references 
therein. Here we just extend them to encompass the class of time varying ph(t) and G(t) models under 
consideration. 
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Here H(t) is the Hubble function, and u(x, t) is the comoving peculiar velocity. We have 
introduced also a perturbation for Gn, equation ( A.7 ), which is mandatory in order to 
have a consistent picture in the model under consideration [45]. Furthermore, in view of 
the corresponding EoS for vacuum and non-relativistic matter, 5p\ = —Sp\ and 5p m = 0. 
Our next task is to insert Eqs. (A.4)— ( ATJty into Eqs. ( A.1| )-( [A~3| ). To this end we first use 
the definition of the gradient with respect to co-moving coordinates, V = V# = a(i)Vf, 
and in this way we can express the result as follows: 



6 m + (1 + «5 m )V -u = 




' U + 2Hu + (u-V)u + -fVd) 



dt a 
AirGa 2 (p m 5 m - 2p A (5 A ) - iira 2 5G (p m - 2p A ) = 



(A.8) 
(A.9) 
(A.10) 



Note that in order to get ( |A.8| ) we have assumed the condition V<5 m = 0, which holds for 
the spherical collapse of a top-hat distribution |33| . Next we take the divergence of the 



Euler equation ( A.9|) while using the following identity 



12 



V-[(u-V)u] = -(V-uf 



(A.ll) 



together with the time derivative of the continuity equation (A.8). Combining all three 
equations, we finally obtain the fully non- linear evolution of the matter density contrast: 



S m + 2HS n 



4 5 2 



4vrG(l + 5 m ) (pmSm - 2p A 5 A ) - 4ir5G(l + 5 m ) (p r , 



3 1 + 8 m 

Inserting now Eq. (|]3|) into ( |A7l2|) , we are left with: 

4 %n 

5 m + 2H5 m - - m 47rG/9 m (l + S m ) 



Sm + 



5G_ 
~G 



. 



2p A ) = . 
(A.12) 

(A.13) 



Changing the independent variable from cosmic time t to the scale factor a through the 
relation dt = aHd a , equation (A.13) can be rewritten 



5" 



— + "77 1 S' m 
a H 



4 5' 2 



3 1 + 5„ 



%^(i + <W 



5G 

m + "77 



(A.14) 



where Cl m (a) is defined in ( p.15 ) and /' = d a f for any /. In the particular case when G is 
constant there are no perturbations of G in (|A.7|) and then equation ( A. 14 ) boils down to 
Eq. (18) in [fjSfl , as it should (see also Eq. (7) in [j33"f ). Indeed, for constant G, the relation 
d3.15[) tells us that Q m {a) = Q m {a)/E 2 (a) = Q m / (a 3 E 2 (a)) . Let us also remark that 
flA.14 ) reduces to 



S m + 



- + — ) 

a H 



30 m (a) 
2a 2 



Sm ~\~ 



6G 
~G 



12 Note that in (A.ll) we have assumed vanishing shear and rotation tensors 
spherical symmetry with a top-hat profile. 



, (A.15) 

owing to the assumed 
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if we keep only the linear terms. The latter is formally identical to equation ( 3.32 ) of 
Sect.g, 

Next, following the prescriptions of fl38| , we compute 5 c (af) (for any scale factor aj, 
in particular for aj = ao = 1) by numerically integrating Eqs. ( [A, 14 ) and ( |A.15| ) in the 
following manner: 



1. First, we run the second order non- linear differential equation ( A.14| ) between Oj and 



df (where ai is a sufficiently small scale factor, which we take as 10 -6 ). Our aim 
is to find the initial value <5 m (oj) for which the collapse takes place at a = aj, i.e. 
such that 8 m (a,f) is very large (formally infinite) at the collapsing time. In practice 
(in order to set the initial conditions), we may assume that the collapse is achieved 
once 5 m is sufficiently large and set e.g. 5 m (af) = 10 . This coincides with the value 
chosen by |38| , although we have checked that our results remain practically the same 
if we take different large values for 6 m (a,f), say 10 5 or 10 9 . As for the initial condition 
on S' m , and since we know that at a = ai this derivative should very small, we can 
take (once more following |38|) 5' m (ai) = 5 • 10~ 5 . Again, any other small number 
(including 0) would yield virtually the same results. 

2. Second, we use the value for 5 m (ai) computed in the first step (together with 5' m (ai) = 
5 • 10~ 5 ) as the initial condition for the linear equation (A. 15). Solving the latter for 



r. 



5m(a>f) we find 5 c {af) by definition of this quantity (the linearly extrapolated density 
threshold above which structures collapse). 

Before proceeding with the above two-step procedure, we have to get rid of 5G on the 
h.s. of ( |OD and flOS|) as follows. Since the growth factor in a pure matter (Einstein-de 



Sitter, i.e. CDM) universe evolves as -DEdS = a i we normalize our growth factor such as 
to get D ~ a at early enough epochs due to the dominance of the non-relativistic matter 
component. Thus we first solve the third-order linear equation ( 3.34] ) between aj and df. 



We take D{ai) = di, D'{ai) = 1 and -D"(a,) = as the initial conditions, and then use 
the solution to construct an interpolation function for SG(a) through ( p.33 ) [in order to 



do that we assume <5G(oj) = 0]. Using this function 5G(a) as an input for Eqs. ( A. 14 ) and 
JA.151 ) we finally proceed as explained in steps 1) and 2) above. The typical behavior of 
5G(a) can be seen in flql , Fig. 3b. 

In the left panel of Fig. || we present the evolution of 5 c (z) in the A^GfCDM model, 
using the different values of v that have been considered throughout this paper. The corre- 
sponding values at the present epoch (effectively defining the collapse time) are indicated 
by S c = <5 C (0). They were used in the number counts analysis of Sect .|7.2|, and are listed 
in Table 1. In the left panel of Fig. |6| we have also included 5 c (z) for the standard model 
ACDM (for which 5 C (0) = 1.675), and also the constant value for the CDM model 13 . 

For comparison, we have also computed S c (z) by neglecting the perturbations in G(a) 
and pa (a) (5G = 5a = 0). This would correspond to the canonical or simplest approach 



(i) mentioned in subsection 7.1 . Although the variable G{a) and Ph{o) affect now the 



dynamics only at the background level, they do influence non-trivially the evolution of 



5 For spherical collapse, the CDM result is known and can be computed exactly: S c — ^ (12 ty) 2 ^ 3 ~ 1.686. 
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Figure 6: The critical overdensity 8 c (z) as a function of the redshift for the different models 
considered in the paper. The CDM and ACDM determinations of S c (z) (with Sl^ = 0.284) are 
represented by the dashed gray and solid black lines respectively. For the A t GtCDM model, we 
study the case where we allow for perturbations in pa and G (left panel), with = 0.283 and 
three different values of v (y = —0.004, black squares; v = —0.001, black crosses; v = 0.002, green 



triangles) - cf. Table 1 of Sect. 7.2. Furthermore, for this model we also consider the situation 
where these perturbations are neglected (right panel). Finally, the dashed red line in the left panel 
corresponds to the best fit value for the A t CDM model (i.e. ^ = 0.284, v = -0.0017). The precise 



values <5 C (0) for the two RG models with different inputs are also collected in Table 1 of Sect. 7.2 



the matter perturbations through the modified Hubble function. Indeed, in this approach, 
we have to solve equations ( [A. 14 ) and ( A. 15 ) for 5G = 0, with H{a) given by ( 3.16j ) - 
where both G and p\ are variable. Alternatively, we may use formally the same approach 
(in which 5G = and G = Gq) if we employ our equations (|5.l| )- (|5.2|) and the 



as in 

non-trivial effective EoS given by (|5.4|) . The corresponding solution for 5 c (z) is displayed 
in the right panel of Fig. 6. We can see that there are some differences with respect to 
the results obtained in the left panel for non-vanishing perturbations of p\ and G, but the 
numerical deviations are not dramatic (at the few per mil level). 

The A t CDM model: 

As before, we want to derive the second-order differential equation governing the non- 
linear evolution of the matter perturbations. In this case, however, the continuity equation 
involves exchange of energy between matter and vacuum: 



dt 



+ • (p m v) 



-PA 



(A.16) 



The Euler and Poisson equations ( |A.2| ) and ( A. 3 ) remain unchanged, although we should 
remark that the Newton's coupling G/v = Go is now strictly constant. As before, we 
introduce comoving coordinates and define the perturbations as in (A. 4) - ( |A.6j ), but 



now (consistently with the approach we took in the paper) we will neglect the potential 

t dens 
and (p 



perturbations in the running cosmological constant density pa [so in (A.4) we set % 

(DT 



m 



Inserting Eqs. (A.4)-([A^) into Eqs. (A.l 



we have: 



S m + (1 + 5m)V • U = -Q(t)6„ 



(A.17) 
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du , 1 -> 

— + 2Hu + (u ■ V)u + -rV^ 

m 



, 
, 



(A.18) 
(A.19) 



where Q(t) was defined in Sect.|3], see equation (4.11). By taking the divergence of the 
Euler equation flA.18 ) [using once more the identity ( A.ll )], the time derivative of the 
continuity equation ( A.17| ) and combining all three equations, we arrive at the fully non- 
linear evolution equation: 



2x2 



5 m + (2H + Q) 5 n 



A5j + 5Q5 m 5 m + Q*8, 
3(1 + 5 m ) 



+ 



0. 



where we also used the background continuity equation ( 
regime (A. 20) reduces to equation ( [4. 10| ) of Sect.|4|: 



2HQ + Q - 4irG p m (l + S m ) 

(A.20) 

. As expected, in the linear 



5 m + {2H + Q)5 m - 4vr G p m - 2HQ - Q 



0. 



(A.21) 



It is convenient to change the independent variable from t to a and perform the numerical 
integration of ( A.20| ) and ( |A.21 ) following the two-step procedure 1) and 2) described in 
detail above. In particular, equation ( A.21] ) can be rewritten in the scale factor variable as 



S'Jn + 



3 H' 
a + ll aH 



Q 



2 



.(«) 



2Q 
H 



H 



Q' 



a 2 



0. 



(A.22) 



where again /' = d a f for any quantity / in this equation. For the particular case of 
time independent vacuum energy (p\ = 0) we have Q = Q' = 0, and then upon using 
O m (a) = Q° m / (a 3 E 2 (a)) we see that QOp boils down to equation (19) of Ref. ||, as 



it should. The solution of ( A.22 ), as a part of the aforementioned two-step procedure, 
provides 5 C as a function of the scale factor or as a function of the redshift z = (1 — a) /a, as 
shown in the left panel of Fig. ^ (see the dashed red line in that figure) . For the numerical 
analysis we used the inputs Q, m = 0.284 and v = —0.0017, as indicated in Table 1. The 
corresponding result for 5 C = 5 c (z = 0) can be directly read off Fig.|(| and is quoted in that 
table. 

From the range of 5 C values obtained in Table 1 for the two models, we see that they 
vary from 1.671 to 1.685, and hence correspond to variations of —2.4 per mil and +5.9 per 
mil, respectively, compared to the ACDM model value (<5 C (0) = 1.675). 
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